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Abstract 

We formulate the functional Bethe ansatz for bosonic (infinite dimensional) representations of 
the Yang-Baxter algebra. The main deviation from the standard approach consists in a half infinite 
Sklyanin lattice made of the eigenvalues of the operator zeros of the Bethe annihilation operator. By 
a separation of variables, functional TQ-equations are obtained for this half infinite lattice. They 
provide valuable information about the spectrum of a given Hamiltonian model. We apply this 
procedure to integrable spin-boson models subject to both twisted and open boundary conditions. 



^f^ i In the case of general twisted and certain open boundary conditions polynomial solutions to these 

> 

TQ-equations are found and we compute the spectrum of both the full transfer matrix and its 

00 " 

^O ■ quasi-classical limit. For generic open boundaries we present a two-parameter family of Bethe 

in 

^y-^ . equations, derived from TQ-equations that are compatible with polynomial solutions for Q. A 

o, 

i^ , connection of these parameters to the boundary fields is still missing. 
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I. INTRODUCTION 

The method of separation of variables is a technique reducing a given multidimensional 
spectral problem to a set of uncoupled one dimensional equations. Although this reduction 
can in principle be applied to any integrable eigenvalue problem, its realization is a math- 
ematically hard problem. In this paper we focus on quantum integrable systems provided 
by the Quantum Inverse Scattering Method (QISM). In this framework, representations of 
a quadratic Yang-Baxter algebra allow for the construction of Hamiltonians from a family 
of commuting operators generated by the transfer matrix. To actually solve the eigenvalue 
problem for these operators Bethe ansatz methods are applied. For physical problems with a 



m 



2|is the method of choice for this step. 
^ require the knowledge of a simple 



f/(l)-symmetry the algebraic Bethe ansatz (ABA) 

This approach and generalizations as presented in y 

known eigenstate, the so-called reference state or pseudo vacuum. Unfortunately, it is the 

identification of this reference state that is severely hampered in the absence of total-number 

conservation. As a consequence, alternative methods are needed for the computation of the 

spectrum that avoid this difficulty. 

Many of the existing alternatives, most notably Baxter's method of commuting transfer 
matrices [a] and Sklyanin's functional Bethe ansatz (FBA) [g, |7|, are based on analytical 
properties of the model due to their construction within the QISM and implicitly encode the 
eigenvalues into solutions to certain functional relations. Various flavours of such functional 
methods have been successfully applied to models where no reference state was known, 
including systems based on non-compact symmetries (e.g. the quantum Toda chain |8| and 
the sinh-Gordon model |9|, [lO|) or spin chains where non- diagonal boundary fields break the 
t/(l)-symmetry underlying the applicability of the ABA 4, lllMl7l|. 

Unlike the ABA, this approach does not rely on the a priori knowledge of a reference 
state. Instead, the the representation of the Yang-Baxter algebra underlying the integrable 
model is dealt with on a functional space isomorphic to the Hilbert space of the model. 
This allows to formulate the many-body eigenvalue problem in such a way that it can be 
separated into equivalent one dimensional equations. This 'separation of variables' 43| then 
leads to the functional (so-called TQ-) equations mentioned above. 

In this work we apply the FBA to models that include a bosonic representation of the 
Yang-Baxter algebra for the case of 3^[su(2)]. Particular emphasise is given to integrable 



spin-boson models with a manifest violation of the U{1) symmetry; but also U{1) symmetric 
model will be considered. 

Interactions between spins and bosonic degrees of freedom are an archetypical problem 
in many areas of physics. They are encountered in the description of systems as diverse as 
impurity atoms in condensed matter on top of a phononic background, dissipative quantum 
systems, and of course all sorts of systems involving the interaction of matter and radiation 
in quantum optics. In particular the class of systems, where atoms or ions are trapped and 



controlled for various purposes [18|, Il9| , has experienced a boost of interest in the context 
of quantum information processing and the simulations of complex physical systems. Many 
of the latter address problems from condensed matter physics, although there exist also 



approaches to study quantum field theories this way (see e.g. Chapters 6&7 of Ref.[l8|]). The 
spin-boson interaction can be decomposed into what is called a rotating part, i.e. aS~^+a^S~ , 
and a corresponding counter-rotating part aS~ + a'^S^. In the presence of only one of these 
terms, say the rotating part, the interaction leaves the U{1) charge S^ + neosons conserved. 
In this case, the Hamiltonian model is block-diagonal with separate blocks for each value of 

nn 

the conserved f/(l) c harge I20l422j. Many integrable extensions of these models have been 
found and analyzed 23|-|26j, where an integration of certain non-linear interaction terms 
has been achieved while keeping the ^(1) symmetry. Hamiltonian models including both 
rotating and counter-rotating terms have been obtained in Ref. [23] using the QISM, and 
by imposing suitable open boundary conditions. 

The article is organized as follows: In the next section we sketch central elements of the 
QISM and give a concise introduction to the functional Bethe ansatz for systems subject to 
quasi-periodic (twisted) and open boundary conditions. In Section Hill we apply the FBA to 
models on an infinite dimensional Hilbert space describing both spin and bosonic degrees of 
freedom. Again, we consider different boundary conditions: it is known that for arbitrary 
quasi-periodic boundary conditions the spectrum can be obtained using the ABA. Here we 
reconsider this case in an FBA approach in Section IIIIAi Open boundary conditions are 
relevant for the models with both rotating and counter-rotating terms in the hamiltonian 
from Ref. 27| . In Section IIIIBI we present in detail the derivation of the TQ-equations 
for this case. The spectral problem for the spin-boson model with both types of boundary 
conditions as encoded in these functional equations is investigated in Section IIV[ In either 
case we consider both the full transfer matrix and its so-called quasi- classical limit. The 



procedure for taking the latter is sketched in the appendix, it extends the family of integrable 
hamiltonians which can be obtained within the approach used for the construction of the 
model. Within the FBA we find a complete solution of the eigenvalue problem in terms of 
a set of algebraic 'Bethe' equations for boundary conditions which can also be treated using 
the ABA. Within the generic set of open boundary conditions leading to both rotating and 
counter-rotating interaction terms in the hamiltonian we propose a three-parameter family 
of Bethe equations, which are obtained from the TQ-equations but using a factorization of 
the quantum determinant that allows for polynomial solutions for the function Q. 

II. FUNCTIONAL BETHE ANSATZ 

The functional Bethe ansatz method was originally formulated as a constructive way to 
realize a separation of variables of a many body system, namely reducing a multidimensional 



problem to a suitable set of one dimensional ones[6j. The method relies on the concept of 
quantum integrability as provided by the Quantum Inverse Scattering Method [l]. It may 
give insight to the exact spectrum for those systems where the ordinary algebraic Bethe 
ansatz fails 8|. In this section we sketch the functional Bethe ansatz for quantum integrable 
systems of interacting spins. 

The basic object of the QISM is the quantum i?-matrix satisfying the Yang-Baxter equa- 
tion 

i?12(A)i?l3(A + A')i?23(A') = i?23(A')i?13(A + A')i?12(A) . (2.1) 

It acts on a tensor product V ®V ®V of a vector space V^ of a given dimension as a function 
of the so-called spectral parameter A. The indices of Rij indicate on which copy of the 
tensor product the i?-matrix acts non-trivially. The i?-matrix can be exploited to provide 
the commutation rules of an associative algebra T (quantum affine algebra) as 

i?i2(A - \')T^^\x)T^^\>:) = r(2)(A')r«(A)i?i2(A - a') , (2.2) 

where T^^-* = T® H and T^^-* = ]1®T and T can be considered as an operator- valued matrix 
of dimension dim V. The vector space V is referred to as the auxiliary space. The algebra 
T is of relevant interest in the theory of integrable quantum systems because each of its 
representations provides a family of commuting operators. From this family a hamiltonian 



model is deduced and the members of the family can then be considered as integrals of the 
motion. 

In this paper, we will consider exclusively the rational solution of the Yang-Baxter equa- 
tion identifying the Yangian affine algebra 3^[s/(2)] 
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A. Quasi-periodic boundary conditions 



Quantum models with periodic boundary conditions are constructed within QISM by 
choosing the representation of T as 



7{X) = CLiX)CL-liX)...C^iX). 



(2.4) 



The matrix Cj is the so called Lax matrix. It is of dimension dim V and has operator- valued 
entries acting non-trivially in the quantum space of site j only. These Lax matrices also 
have to fulfill a Yang-Baxter equation 



'{!)/ 



'(2)/w 



(^)f\'\rW, 



Ri2{x - x')c)'>{x)cy>{x') = cy>ix')cy>{x)Ru{x - x' 



(2.5) 



with the additional requirement of ultra locality Cj (A),£j (A) =0 for j ^ I. Quasi- 
periodic boundary conditions can be used as the simplest way to introduce boundary terms 
to the final hamiltonian. In the realm of the QISM this can be done slightly modifying T as 



-(L) 



TZiitW = KT{X) = 



A{X) B{X) 
CiX) D{X) 



(2.6) 



We comment that for the rational i?- matrix from Eq. (12. 3p the relation (12. 2p is satisfied by any 
A independent C-number matrix K of dimension dim V in (12. 6p . since [i?i2(A), K'^^^K^'^^] = 0. 
The generating functional for the hamiltonian, and as well for the integrals of the motion, 
is the transfer matrix ttwist(A) = try7twist(A) with the trace taken over the auxiliary space. 
The FBA method allows to construct separation of variables for the spectral problem 



ttwist(A)|^)=A(A)|7A) 



(2.7) 



The starting point of the procedure consists in looking at the operator valued zeros of the 
'lowering operator' C(A) enjoying the property 

[C(A),C(/i)] = 0, VA,/i (2.8) 

(equivalently, one can choose to consider operator zeros for the 'raising operator' B{X) with 
[-B(A), B{^)] = , VA, /i). For K21 7^ 0, the meaning of the roots of the operator C(A) can 
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L 



be specified by expressing the latter as J6| 

C(A) = K,, J] (A - x„) (2.9) 

n=l 

where the operators x„ can be simultaneously diagonalized because of the vanishing commu- 
tator [xn, Xm] = descending from the basic commutation relation (12. 8p . In turn we observe 
that C(A) can be diagonalized as it is indeed a polynomial operator of order L in the spec- 
tral parameter with coefficients that are symmetric functions of the roots Xi,X2, ■ ■ ■ ,Xl- 
Therefore the operators A{X) and D{X) are not diagonal in this basis. We define 

A{\ = Xn):=J2xlAp = X- 

" (2.10) 

p 
where an operator ordering is established by placing x to the very left in each term. We define 
Sym[xi, . . . , xl] as the set of symmetric functions of arguments xi, . . . ,xl- The operators 
X^ act on elements of Sym[xi, . . . , xl] as 

X± Sym[xi, ...,xl]= e±^^/^"" Sym[(xi, . . . , xl)] = Sym[xi, . . . , x„ ± r^, . . . , xl]X^ (2.11) 

suggesting them to be considered conjugated to the operators x„. In fact, the commutation 
rules are |6| 

[A^,Xn\ = ±TlA^Omn 

On a generic (not necessarily symmetric) function f{xi, . . . , xl) the operators X^ act as 

XnfiXl, ...,Xl) = A^{Xn)fiXi, . . . , X„ ± ?], . . . , Xi) (2.13) 

where A^ provide a factorization of the so called quantum determinant detg of the mon- 
odromy matrix Ttwist- The quantum determinant can be expressed as detg(7twist) = 



A{\ + ri/2)D{\ - ri/2) - B{\ + r//2)C(A - ri/2) yielding directly 

X^X^ = A±(x„)A^(x„ ±r]) = det,(7;„ist(x„ ± vf^)) ■ (2.14) 

Evaluating by substitution from the left the spectral equation (12 .yp at A = x„ and sand- 
wiching between left and right eigenvectors i{x\ and \x)r of xi, . . . , £„ gives 

A{xn)i'{xi, ...,xl) = S+A+(x„)7/'(xi, . . . , Xn+V, ■■■, xl)+E:" A"(x„)V^(xi, . . . , Xn-r], ■■■,xl) 

(2.15) 



where '?/'(xi, . . . , xl) = i{x\il>) and E^ = (tr i^i ^/(tr K)"^ — 4 det K)/2. The final separation 
of variables is achieved by the ansatz ip{xi, . . . , xl) = Ylj Qji^j) leading to 

A{Xn)QniXn) = S+ A+(x„)(5„(x„ + T]) + S" A" (x„)Q„(x„ - T]) . (2.16) 

These equations are L one dimensional finite-difference equations for a;„ G with & being 
the 'lattice' provided by the eigenvalues of Xi, . . . , x„. 

B. Open boundary conditions 

The FBA method has been generalized recently to integrable models with open bound- 



aries in Ref. [iTj]. In the seminal paper [28|, |29| Sklyanin demonstrated how to enlarge the 



class of integrable models obtainable from QISM by defining the so-called double-row transfer 
matrix describing a closed system interacting with a boundary. The core of the construction 
is the set of reflection algebras 

Ri2{\ - \')K^l\\)R2i{\ + A')i^i'^(A') = iri')(A')i?i2(A + A')i^i'^(A)i?2i(A - A') , 
i?2i(-A + A')i^l'^*(A)i?i2(-A - A' - 27^)Kf\\') (2.17) 

= irf *(A0i?2i(-A - A' - 2r^)Kf{X)R,2{~X + A'), 

where K{X) parameterizes the boundary conditions, and K^ is the transpose of K. The 
involved i?- matrix is again a solution of the Yang-Baxter equation (12. ip . Additionall y it 



fu 



fills the conditions of unitarity, parity, time reversal invariance, and crossing symmetry 28 



29| . It can be demonstrated that the following objects are representations of the reflection 



algebras f ETfp : 

r(+)(A) = ir+(A) , T^-\\) = T{\)K_{\)T-\-\). (2.18) 



Within this framework the generating functional for commuting integrals of the motion in 



Eg. ( 12.71) is the following operator (double row transfer matrix) 



i?e.W=tTvT^^\X)T^-\X)=tTy [K^{X + vmX)K4X)T-\-X)] . (2.19) 

As customary, we define W(A) = detq[T(— A — r]/2)]T^^\X) with its matrix representation 
on the auxiliary space 

W(A) = I ^^^ ^^^^ I . (2.20) 

IC(A) ViX) 



We particularly note that [B{X),B{fj,) 



(for further relevant commutation relations be- 



tween the operators in U{X) see e.g. 28|). In case of the rational i?-matrix from Eq. 02.3p 
the general non-diagonal C-number representations of the reflection algebras are the K- 
matrices [SoJ K+{X) = \K{X + r], +) and i^-(A) = K{X, -) with 



^(A,±) ,i 



1 / A + e^ 2fi:V*A 

^"^ \2«±e-^*A -X + e 

(2.21) 
1 X sinh (3^ + a^ cosh P^ Ae^ 



a±cosh/3± I ^g-e± -Asinh/3± + a±cosh/3± 



where (see Ref. HJ]) a^cosh/3^ = ^^ and sinh/3^ = y^- ^^ the latter parametrization a 



diagonal boundary corresponds to the limit /?''=—;■ oo. 

For the double-row transfer matrix (I2.19P the FBA method proceeds through similar 
steps as discussed above for the quasi-periodic case but applied to the matrix W(A) instead 
of %l,iltiX). B{X) in terms of its operator zeros can then be expressed as 

The property [a;f,x^] = 0, arising from the commutation relations of the operators in 
U{X), can be assumed as emerging from [a;;,Xm] = 0. We sometimes write formally xi = 
diagja;/', x^} in terms of the eigenvalues x^ of the non-hermitean operators xi. The operators 
V{X) = 2XV{X)-r]A{X) andC{X) = (2A -r/)C( A) give rise to the shift operators X+ = A{xn), 
and X~ = 'D{xn)- The action on generic functions is given in Eq. (l2.13p as the operator 
valued zeros provide the same algebra (I2.12p and provide the factorization of the quantum 



determinant detqU{\) = ^(A + t//2)P(A - r//2) - B{\ + r//2)C(A - ri/2): 

X-X+ = A-(x„)A+(x„ -r]) = det.Uixn - r//2) 

X+X- = A-{xn + r/)A+(x„) = det,U{xn + v/2) 

where 



(2.23) 



xf) 






(2.24) 



AL) 



The spectral equation for topen(A) from Eq. (l2.19p reads as 

KMij{x,, ...,xr.)= ^-^^^^^^^^^p^±^^^^^A^x^)ij{x„ ...,x^ + v,...,xl) (2.25) 

x„-a+ + r//2 

H T ^^ A (x„)V'(Xi, . . . , Z„ - 7], . . . , Xi) . 



■-n'- 



The ansatz '?/'(xi, . . . ,xl) = Ylj Qji^j) again leads to an uncoupled set of one dimensional 
finite difference equations 

A/ x^ / X {Xn + V/'2){Xn + a+ - v/2) ^ 
A{Xn)Qn{Xn) = ' ^ A^{Xn)Qn{Xn + V) 

+ -^^-^ + A'{Xn)Qn{Xn - Tj) . 

When we extend the validity of these TQ-equations to complex x, only a single function 
Q{x) remains. This will be our working hypothesis. 



III. FUNCTIONAL BETHE ANSATZ FOR SPINS INTERACTING WITH A SIN- 
GLE BOSONIC MODE 

In this section we will carry out the functional Bethe ansatz for spin-boson hamiltonians. 
This corresponds to generalizing the approach to an infinite dimensional quantum space. In 
the framework of the QISM, integrable models for interacting bosons (e.g. photons) and 
spins (e.g. two-level atoms) have been constructed from the algebraic structure induced by 



m, 



27| using the boson and spin 



the i?-matrix (12. 3p for the rational six-vertex model 12 
Lax operators 

, A — rizi — ria^a Ba^ \ /A — nzn + nS^ rtS' 

aA)= I ,1; Cs{X)= \ ' ' I (3.1) 

la -^ / \ r]S+ X-tjZq- r]S' 



where 2:05 ^i cire inhomogeneities for the spin and the boson, respectively. The quantum 
determinants are detg (£5) (A) = —^{X—{zi — ^)ri) aB.ddetq{Cs){X) = {X—rizo—ri){X—T]Zo+ri); 
the bulk monodromy matrix is defined as 

r(A) = CbiX)Cs{X) . (3.2) 



A. Twisted boundary conditions 



E 



The FBA for models with twisted boundary conditions has been introduced in Ref. iGl 

n n 

and applications include e.g. the Gaudin model [7[ and the Toda chain [8|. Here we apply 
this technique to models including a single bosonic degree of freedom. An extension to more 
spins and/or bosons is straight forward. 

Spin-boson models derived from the XXX i?-matrix defined in Eq. fl2.3p and with twisted 
boundary conditions are known to lead to models without counter-rotating terms [27!] • This is 
a consequence of the observation that every boundary twist matrix can be brought into upper 
triangular form by means of .oca. gauge transformations g g where the transformation 
matrices are elements of the symmetry group of R [SJJ. Although diagonalizable via the 
ABA, the application of the FBA to these simple models is still interesting for two reasons: 
at first, we can compare with the results obtained from the usual algebraic Bethe ansatz 
approach, and second we will demonstrate explicitly how the FBA machinery works when 
bosonic degrees of freedom are included. 

We start from the monodromy matrix with general twist matrix {K)jk 

T^W = I ^'' ^'' I UX)C,{X) (3.3) 

\ K21 K22 J 

and find that B has no term ~ A^. Instead, such a term is contained in C, and we will 
perform the FBA and separation of variables method for the operator C instead of B (this 
"asymmetry" is a consequence of the peculiar form of the bosonic Lax-operator). 

It turns out to be convenient to consider linear combinations of B and C rather than 
C directly. This can be realized by suitable similarity transformations of the monodromy 
matrix which do not affect the resulting transfer matrix. We summarize this procedure in 
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the monodromy matrix 



a (5 \ h- b\ A B 

A(A)£,(A) 
1 c* / \ c d \ C D 



(3.4) 



In complete analogy to Ref. l32| we define the operator zeros Xj of C as 



N 



N 



c{x) = n(A 



Xj) 



j=o 



(3.5) 



After subsequently applying the shift a — )■ a + c to the bosonic operators for convenience, 
we find 



ci = -r] [Sz + cS - n] + zo + zi 

C2 = —rf \_h{Sz + cS~) + a){2cSz + c^S^ — S^)\ — rfZiCi + {zi — ZQ)h — zqZi . 
A common basis of right and left eigenstates of both operators is spanned by 



(3.6) 



2y/m + l 



\+,m) = I i)\m) + - , 

' ' ^ '^^' ^ 2{m + zi - zo) + 1 

\-,m) = {\t) + c\i))\m) 

(+,m| = -(c(t|-ai)(m| 

2./m 



(I t) + cU))|m + l) for m > 



(3.7) 



[^2,: 



(t|(m| + 



2(m + zi - zq) - 1 



'c(t I - (i l)(^-l| for m > 



±,m| 



and the corresponding eigenvalues defined by Cj\±,m) = c- '"^{±,171) are 

Ci'"" = ri{m + Zo + Zi ± ^) ; €2'"" = r]'^{zo± \){m + Zi) . 



(3.8) 



The spin and boson operator zeros Xg and Xb are uniquely determined from the spectral 
decomposition as 

1 



Xs = -r] 



2a^- 



Xb = V 



2{n + zi-zo) + l 
1 



+ ^1 + 2a'^ 



n + zi 



2{n + zi- zq) + 1 



{2cSz + c^S- - S+) + cS- + Sz~Zq 
(2cSz + c^S- - S+ 



(3.9) 



and their eigenvalues are 



Xs,± = ri{zo±\) ; Xb^rn = ri{f^ + Zi) . 



(3.10) 



In the limit of a diagonal boundary twist, pairs of the right and left eigenvectors of Xg and 
Xb coincide and the procedure breaks down. 
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Having the spectra of the operator zeros of C, we can exphcitly write down a factorization 
of the quantum determinant detq(T'' '*''•* (A)) = A+(A— 2)A^(A+|) that leads to TQ-equations 
after a separation of variables 

Q{xs,Xb) = Q{xs)Q{xb) . (3.11) 

We find 

A+(A) = -^(A-r^(zo + |)) 

A{x)Q{x) = E:+A+(a;)(5(a; + r]) + E:"A"(x)(5(x - r]) 
with S± from Eq. fl216|l . 

B. Open boundary conditions 

In order to create counter rotating spin-boson hamiltonians we have to resort to open 



boundary conditions |27|, |28| . In this case the FBA can be applied as well [17| as described 
above. 

We consider the following transfer matrix 

^Sn(A) =trir+(A- V2)W(A) (3.13) 

with U{X) defined as 

U{X) = r(A - v/^)K-iX - V2)a^T(-A - V2)^^ = f ^^^^ ^^^^ I • (3.14) 

\CiX) V{X)) 

The boundary matrices are chosen as K_{X) = K{X,—) and K+{X) = ^K{X + ri,+) in terms 
of the ii^- matrix given in Eq. fl2.21j) . 

As i3 is a polynomial in the spectral parameter of maximal degree we directly factorize 
it in terms of its operator valued zeros. In Ref. Il7|, Eq. (4.4) it was shown that 



B(A) = _ (2A-.|)_tanhr g_^_(,) ,3 j^, 

where ^symm(A) = B^X^^ + i?2A^ + Bq is an even function of A. An expansion of -Bgymm in its 



operator valued zeros x is 



sinh/3 — sinh(6' ^ ^^ — /? ;.,2 -2\/\2 -2\ 



12 



We anticipate that the C-number zero A = 77/2 gives rise to an additional constraint on the 
eigenvalue of topen(A) but will focus on the operator- valued zeros first. 

Using the symmetric polynomials 61 = — 52/-B4 and 62 = -B0/-B4 we obtain 

bi=xl + xl ; b2 = xl ■ xl (3.17) 



The explicit coefficients 61 and 62 have been computed with FORM 33|. 

The bosonic Lax operator has entries proportional to a single a^. As a^ does not have 
any right-eigenstates but a does we will consider 6{ and b^ in order to work with kets for 
convenience. Expressed as a 2 x 2 matrix in spin space the operator bi reads 

with bosonic operator valued entries 



S++^ = -a^^e-^®- + a^^l^ {z, - ^ + n) + z^ + [z, + n + i)^ 



7] /« V 2 

b^r^ = -2a^ 
V 

SS") = -a^e-^®-^ + ^^^ {,, + 2-^+n)a+ {zo + if + {z, + n + l: 

Tj^ KT] ^ ' 2 



where n = a^a is the bosonic number operator and ^ = $,/r]. The wave function is written 
correspondingly as a two-component vector 

00 

1^) = ^'^na\na) = ^(V^nth)® lt)+V^n;|^)® 14)) 
na n=0 

Considering the eigenvalue problem (b\ — El) \-ip) = 0, where E is the eigenvalue of b\, 
orthogonality leads to two intertwined recurrence relations for the coefficients ijjna- Defining 
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E = E jrf, (5s = 2:0 + |, and 5f, = -^i + | this implies 

= -Vn + 2Vn + le-2^>„+24 + Vri + 1^ U<5b - ^ + n - ^ j ^n+it + (3-20) 

("(5. - i)' + (5, + n)2 - e\ ^„t - 2v/;7TTe-2^>„+u + ^-^{5s- ^ ^„; 
= -2yjn + lz^„+i^ - Vn + 2(n + l)e-2^ i)n+2i + Vn + 1 x (3.21) 

Ub + ^ + n + n ^IJn+n +\{5s+ 2)' + {^b + nf - e\ tPni 

Considering the large n regime (n large as compared to the corresponding eigenvalues) 
each coefficient satisfies a F-function like functional relation. Hence, in order to obtain 
normalizable states the recurrence relation must terminate at some finite boson number m, 
and the eigenvalue can be read off directly from the coefficient of the highest boson number 
state \m) leading to 

[{5, + mf + (6, - lY - E) i,„,, + [^{5. - 0) i,^\ 
[{h + m)2 + (5, + \f - Ej ^„4 J 

There are two possibilities for satisfying this set of equations 

El' = i5h + mf + {S, - ^f ■ iJ^i = ■ z^„t ^ arbitrary (3.23) 

&,' = {6, + mr + {6. + lr ; ^ = f^i^s ~ ■ (3-24) 

2 Ipml 2K,ds 

The corresponding eigenstates can then be calculated by carrying out explicitly the recursion 
setting ipm+i,. = and ipm,. as stated just above. 

An analogous calculation for the operator 62 (with E = E/r]'^) leads to 

El' = {Sh + mf ■ {5s - ]^f ; V'mi = ; V-^t 7^ arbitrary (3.25) 

The eigenstates again result from recurrence relations. 

As hi and 62 commute due to [i3(A), i3(/i)] = they share a common system of eigen- 
vectors. The eigenvalues of hi (62) are only degenerate for a finite number of states if the 
inhomogeneities zq, ^1 (resp. 5s/b) are chosen carefully. Typically, this degeneracy is lifted by 



the 62 (hi) operator. [4^ But e.g. the choice zq = —1/2, i.e. 6s = 0, results in a massive de- 



generacy of both hi and 62, and the separation of variables cannot be carried out in a straight 
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forward way 45| . It is a straight forward calculation that the above right-eigenstates of b\ 
are also right-eigenstates of bl- 

With the eigenvalues at hand it is possible to write the operator valued zeros of B as 
a matrix acting as multiplication operators on the common eigenbasis of b\ and b^ from 
Eg. ( 13.1 7]) . For bosonic quantum number m we find 

As the x^ operators can be simultaneously diagonalized the operator zeros are 

h = v{ \ ; x, = v{ ' 1 . (3.28) 

This also fixes the sets (& for the lattice of the TQ-equations 

&b = {SbV, (Sb + l)v, (Sb + 2)ri, ...}=: X,- + No ; (3.29) 

^s = MSs-l),viSs + l)}=--{x;,xt}. (3.30) 

As discussed above, the operator-valued zeros can be seen as 'coordinates' whose 'conju- 
gated momenta' are the shift operators 

(3.31) 



\=Xi 



v{X)=J2^p,^x] 



where Ap and Vp = 2XV — riA denote operator- valued expansion coefficients of the operators 
of the algebra U and j G {s, b}. These 'conjugated momenta' are representations of an alge- 
bra analog to f l2.12p . On arbitrary functions they act as X^f{xs, x^) = A'^{xs)f{xs ±r/, Xb), 
X^f{xs,Xb) = A^(xb)f{xs,Xb ± 1]) and induce a factorization of the quantum determinant 
of W: A'^jx — f] /2) A^ ( X + v /2) = detqU{x) for all x G 0. The following factorization meets 



these demands (compare [iTj) 



A-(A) = L^llH^^x - r^{6., - 1))(A + ^{6, + b)(— (A - 5,r/)) 

" _ ^ ^ "^ (3.32) 

A+(A) = (2A - r^)^±^^/^^^(A - vi^s + ^))(A + vi^s - ^))(^(A + 5bv)) ■ 

a 2 2 7] 

Note that A^*" vanish on the appropriate boundaries of the sets (3s^b especially the set &b is 

only bounded from below and hence only A^ {SbV) is required to vanish. This is a consequence 
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of the fact that we have an infinite dimensional representation of the algebra fl2.12p . Given 
the lattice & and the above factorization of det^W, we arrive at the TQ-equations 

Aix,)Qix,) =^ 'IJ ^^A ixj)Qix, - T]) 

+ (3-33) 

with an unknown function Q{x). The allowed arguments are Xj G &j on the grid indentifying 

A^±(x) = A^(xj) for J G {s,6}. 

IV. THE Q-FUNCTION AND THE SPECTRAL PROBLEM 

In this section we investigate the possible information on the spectrum of the transfer 
matrix given by the TQ-equations arising from the FBA method. We remark, in this con- 
text, that the function Q in its own does not carry any immediate physical information. 
Instead, the values of Q on the lattice might be the key to the eigenfunctions by means of 
a yet to be found isomorphism. The basic working hypothesis consists in the identification 
of the latter equation with the Baxter TQ-equations, namely by extrapolating the validity 
of the TQ-equations to outside in interpreting them as a functional relation. We comment 
that the FBA method can be viewed at as being complementary to the ordinary algebraic 
Bethe ansatz for the spectral analysis (especially for those problems with no straightfor- 
ward pseudo-vacuum states) in the sense that the method searches for a basis in which the 
'lowering operator' of the Yang-Baxter algebra is diagonal. 

It is worth noticing that the TQ-equations derived on the lattice do not have a unique 
solution in the space of continuous (or meromorphic) functions, in particular not as far as 
the function Q is concerned. Whereas the eigenvalue A of the transfer matrix is known to 
be a polynomial with its degree given by the transfer matrix, we do not have such a priori 
knowledge about the function Q: additional knowledge about the latter is needed as an 
input. 

In those cases where a reference state can be found (possibly after gauge transformations) 
and the algebraic Bethe ansatz can hence be performed, polynomial solutions for Q exist 
and lead to the known Bethe equations. But even in these cases, no general direct link 
between Q and the eigenfunctions is established in the literature. A central open question 
concerns necessary and sufficient conditions for a polynomial representative Q to exist. We 
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will not give an answer to this question, but will discuss specific cases, which hopefully will 
shed some light onto this important problem. 

A. Quasi-periodic boundary conditions 

We start with the TQ-equations from ( 13.121) . We know from previous studies [27| that 
quasi-periodic boundary conditions always lead to algebraic Bethe ansatz solvable spin- 
boson models corresponding to slightly generalized Jaynes-Cummings hamiltonian (without 
counter-rotating terms). According to the general scheme discussed above, we note that 
the continuous limit of the FBA TQ-equations coincides with the Baxter equations found 
previously [12] • Therefore, each class of solutions for the TQ-equations will have a polynomial 
representative: Q{X) = na=i('^ ~ ^a) leading to the Bethe equations 

vKf^{X,3 - yzp - f )(Aff - r]Zi) _ X^-X^ + ri 



h det K-{\p-r^z^~l) J-^-L A^ - A, 



and eigenvalues 



a a 

(4.2) 



Both Bethe equations and eigenvalue equation agree with those in Refs. 12|, |25| 



1. Quasi- classical expansion 



The quasi-c 
formed in Ref. 



assical expansion of the twisted spin-boson transfer matrix has been per- 



12 



usmg 




(4.3) 
leading to an integrable generalization of the Tavis-Cummings model; X, U, V are real pa- 



12 



have been set to 



rameters; U, V have the same sign. The parameters Y and A in Ref. 
zero here. By inspection of the transfer matrix (13.31) with the parameterization (14. 3 p the 
two lowest order contribution in powers of 77 are proportional to the identity. Therefore, the 
first non trivial term in the expansion (resulting oc rj) of the TQ-equations (see IA3p leads to 
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the following Betlie equations of the Gaudin type 



2+ 2VViV-U) '•" ^"^A^^ ^^^ 

after the above polynomial ansatz for Q. These Bethe equations themselves determine a 
second order differential equation for Q: 

XQ" - (1 + ^^^lySj^^ ^ + 2A^)Q'(A) + (MA - C)Q(A) = (4.5) 

where C = -^ = -E.A-^ 



B. Open boundary conditions 

We consider the TQ-equations f l3.33p on the lattice (5. The vanishing of A''^ on the 
boundaries of (S^ = {xj = ri{6s — ^),xf = vi^s + |)} yields the linear system of equations 

./ +.^, +x (x+ +?7/2)(x+ + «+ - ?7/2) ^ _ , ,,^, _, 
A{xt)Q{xt) ='-^ '; I ^ ^^A (x+)g(x, ) 

^ ; (4.6) 

r~ — fy'^ -\- n/2 

Hx:)Q{x;) = - - ^^_J^^ A+(x;)g(x+) 

and A{xf) are obtained from the condition that the determinant of the coefficient matrix 
of the linear system vanishes. The vanishing of A~ at the 'lower' boundary 6br] = x^ of 
^b = {^bV, i^b + 1)^, ■■■,{Sb + n)T], ...} = {xl, xl...,x^,...} resuhs in 

H4)Q(4) - '"° ^ "^"gj;"' " "^'' a-(x;)Q(.,-') (4.7) 

and for the other lattice points we get {n > 0) 



X^ + V2)(X^ + «+ - V2) , _, „,^. (n+l). 



Aix^)Q{x^) = '" ' " "^J^^^ ^A-(xng(x, 



"6 



(4.8) 



Here, A(x^) are formally obtained from a continuant of a half-infinite matrix being zero. (A 
continuant is a determinant of a tridiagonal matrix: see e.g. [34] )• 

The expansion of the eigenvalue A in powers of the spectral parameter A involves 4 
coefficients. As the asymptotics is known, 3 equations are needed to fully determine 
A. Besides the conditions of vanishing 2x2 determinant and half- infinite continuant 
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we have a third equation from the C-number zero 77/2 of the B operator. At this point 
the eigenvalue is the quantum determinant detg T of the periodic monodromy matrix: 
A{ri/2) = detgT{—r]/2). This can easily be seen by considering the unshifted and un- 
sealed transfer matrix t(A) = tT[K^{X)T{X)K^{X)T^^{—X)]. Inserting A = results in 
t(0) = tr[ir+(0)r(0)]ir"^(0)] = tri^+(0) = 1. The transfer matrices are related via 
t(A) = detgr(-A)t(A-r//2). 

In the spirit of the scheme discussed in the beginning of the present section we now 
interpret the TQ-equations f l3.33p as a functional relation which holds for general complex 
spectral parameter: 

A(A)Q(A) = A-(A)g(A - r/) + A+(A)Q(A + r/) . (4.9) 

One finds that the coefficients 

A (A) = ^^ A (a:,) 

- ^^ irA + ?UA + «^-?UA + a--?U 



2r]a+a- A V 2/ V 2/ V 2 

A + ^{Ss + \)) (a - r/(5. - ^)] (A - r/55) 



(4.10) 



A+(A) = A-(-A) 

behave asymptotically as A^ for large values of the spectral parameter. From the definition 
(I3.13P of the transfer matrix the eigenvalues for non- diagonal boundary conditions grow 
asymptotically as 

In the Bethe ansatz solvable case (diagonal or triangular boundary matrices) the coefficient 
of A^ in (14. lip vanishes and the asymptotic behaviour is oc A^. Since in A^ the leading order 
disappears as well, the TQ-equations (14. 9 p can by solved using an even polynomial ansatz 



27. 



for Q{X) in agreement with the Bethe ansatz analysis in Ref. 

In the non-diagonal case the mismatch in the asymptotics for a chosen factorization 
of the quantum determinant must be compensated by the Q-function. Hence, Q will be 
transcendental. On the other hand, we might want to insist on polynomial solutions for Q, 
but then the factorization must be modified accordingly [46|. 

In order to allow for a polynomial form of Q, our intention is to absorb a non-polynomial 
part in a function J-", which in turn modifies the factorization of the quantum determinant. 
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We hence define 

Q(A)^^(-)Q(A) (4.12) 

V 



such that Q{X) is polynomial in A 47|. To this end we assume that J-" satisfies the functional 



relation ^-"(2; + 1) = p{z)J^{z), where p{z) is some rational function. This leads to the foUow- 

Theorem 3.4). 



— "^ j^(— ±1) — 

ing modified factorization of the quantum determinant A='=(A) — > A='=(A]_^= % ^''^(A) 

■'(7;) 



(representing an 'algebra isomorphisms' as mentioned by Sklyanin in Ref. |6 
On the other hand, this redefinition of Q can be seen as a consequence of a change in 
normalization for the right- and left eigenstates of B, which has no physical significance. 
Therefore we can redefine Q such that the asymptotics of the TQ-equations is fixed already 
by A(A) and A=''(A), hence allowing for polynomial solutions Q. The ansatz 04. 121) leads to 
the modified TQ-equations 

u - - 1 ~ 

A{u)Q{u) = A+iu)p{-)Q{u + r^) + A-iu)—^—Q{u - r/) (4.13) 

= A+{u)Q{u + r]) + A-{u)Q{u - r]) . (4.14) 

There are two possibilities for matching the asymptotics. The first is to reach deg[A"('u)] = 
deg[A+('u)] + 2 ('deg' is the polynomial degree). Then p{u) x u^^, and the ansatz p{u) = 
{pooUTj + x + Poo)'^ leads to 



J^{z) = Q{z)r{-z - C - 1) exp{-az) 

X 



exp a = -pooV ; ( 



(4.15) 



PooV 

where Q{z) can be an arbitrary function with period 1. It does not affect the values on &. 

In the above TQ-equations (14.91) we need p^ = « cosh /3+ cosh g- ^^^ ^ must be chosen 
suitably; one condition to be met would be that no spurious upper bound for the bosonic 
spectrum is created. 

The other possibility is deg[A^('u)] = deg[A+('u)] — 2. In this case, p{u) x u, and the 
ansatz p{u) = {poowq + x) yields 

7{z) = Q{z)r{z + C) exp{az) 

X (4-16) 

exp a = PooV ; C = • 

PooV 

— — u 

The asymptotics of the TQ-equations (14. 9 p is then fixed by the choice poo = — g — ^ 



l3-y cosh /3+ cosh /3~ ' 
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A formal polynomial ansatz Q{\) = na=i('^ ^ Aa)(A + Aq,) leads to the Bethe equations 

^"(A/j) = - TT ^/^ ^ ^" + ^ ^/^ + ^" + ^ C4 17) 

A+(-A,M^)p(^ - 1) i-/^ A, - A. - r^ A, + A. - ry 

in both cases. Inspection of the transfer matrix suggests that the full degree A^ is present in 
the element corresponding to V. Assuming that this corresponds to the proper factorization 
for non-diagonal boundaries, we insert the explicit linear form for p{u) (the second case 
above). The Bethe equations derived from this assumption are 

(A^ + ^)i\(s + a+- i)(A^ + a-- 2)(A;3 - VJ^^ - |))(A^ + riJSs + |))(A^ - riS,) 
(A^ - 2)(A^ - a+ + 2)(A^ - a- + 2)(A^ + ry(<5, - i))(A^ - r^(<5, + i))(A^ + r^S,) 

2 (X , A NM , (/■ TN N TT A/3- A„ + r7A;3 + A« + 77 
= -Poo(A/3 + CV){^P + (C - 1)^) 11 T 7 -^-—\ ■ 

It is evident that the parameter C, plays a crucial role. The case of diagonal boundaries 
is obtained in the limit p^o -^ and C — )■ cxd such that p^oC ~^ x/V- The known Bethe 



equations for the diagonal case [27|, 



28| are obtained if additionally x "^" 1 in this limit. It is 
easy to check that the functional relation for the function J-" reads J^{z + 1) = J^{z) in this 
limit, which is satisfied by constant J-", and hence no transformation of the TQ-equations is 
induced. In case of linear p{z), the simple zero zq = \ — x/iPoo) oip{z/rj — 1) must coincide 
with one of the zeros oi A^{z)Q{z — ri) [48|]. In the TQ-equations this leads to the constraint 



- PoovH-x/Poo)Q{-x/Poo) = A {-x/Poo)Q{-x/Poo - v) ■ (4-i9) 

The correct x "will eventually be encoded in the boundary matrices; it depends on the 
parameters Poo and (. This is also seen from parameter counting: besides the four eigen- 
values, there remain two further parameters. Therefore, the three parameters p^o, (, and x 
can not be independent. 

We emphasize that the factorization proposed above is only one possible choice and 
might not be the correct one. In order to determine the proper factorization, a deeper 
understanding of the FBA is necessary. We leave this open for future research. 
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1. Quasi- classical expansion: diagonal K 

We first consider the TQ-equations f l4.9p for tlie Betlie ansatz solvable case corresponding 
to diagonal boundary matrices fl2.2ip . i.e.: 



K{\) = -r . (4.20) 




For this choice the leading terms of the quasi-classical expansion of the transfer matrix (I3.13P 
in the parameter 77 are 

r(o) = 0, r« = (e+ + r)A^ ^^-^l) 

r(2) = (25, -2n-zi-l)\^+ {i+C{2S, - z{) + 2030)8+) \^ . 
The eigenvalues of r'-^^ in the sector with n — S^ = {2k — l)/2 are 

A^(A) = -(2A: + .0A^-e+^(^l + l)A^ ^ ^ 

(4.22) 

Ag(A) = - (2A: + ^0 A^ - e+r (^1 - 1) A^ 

(2) 

The fc = subspace of the system is one dimensional, the corresponding eigenvalue is A^ a(A). 
At second order in rj (14. 9 p turns into a first order differential equation for the eigenfunctions 

Q{z),z = X': 

z{z + tn Q'Uz) - {kz + en QkA^) = o , ^ ^ 

(4.23) 

{z+enQ'u^)-kQk,b{z) = o. 

Up to normalization these equations are solved by 

QUz)^{z+er)', , , 

(4.24) 
Qfc,6(^)ocz(2 + e+r)'"' for A; 7^0. 

As mentioned above, the Q-functions are even polynomials in the case of diagonal boundary 
conditions - this property still holds in the quasi-classical limit. 

Alternatively, we may solve the TQ-equations without using our knowledge of the eigen- 
values: with the ansatz A*^^)(A) = oaX^ + ^^(^2^ + (^^^~)^Oo we obtain the following 
differential equation for Q{z) (see Eq. flA3p ): 

2z^[z^ec)Q\z) ^^_^^^ 

+ ((a4 + zi) z^ + tC (a2 + ^1 - 1) ;2 + {tCfa^) Q{z) = . 
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The general solution to this equation is 



Q{Z) OC ^(l-.+«0-.i)/2 (^ + ^+^-) (-l-a4+a,-ao)/2 ^^^ (^^\ . (4.26) 

Requiring Q(2;) to be analytic with at most a simple zero at 2; = one obtains immediately 
ao = and the solutions given above is reproduced. 



2. Quasi- classical expansion: non-diagonal K 

Generically, non-diagonal boundary matrices lead to non-hermitean transfer matrices. 
Within the quasi-classical approach, however, it is possible to construct hermitean hamil- 
tonians for the spin-boson model by fine tuning the dependence of the system parameters 
on the 'quantum parameter' rj [27j: rescaling zq — > zq/ii] in Eq. (13. 1|) and parametrizing the 



boundary matrices as 



, e± + A A/i± , 
K{\) = p / (4.27) 



the following choice of parameters 

B B^ (4-28) 

/i+ = r7^(/ir + z/r), ^^+ = , e+ = -^ + e+ + r^^+. 

leads to the following self-adjoint hamiltonian 

H = Qon + As,5, + iA,,(5+ + S') + g{S+a^ + S^a) + 2a{a + a^) (4.29) 

where the coupling constants are obtained as 

n, = 2{zl - A^), A3, = 2(A2 - p\i- - zi)), 

P (4.30) 

Asx = -2/32^0^^1 , g = 2/32;o, a = -u^ (A^ - z^) 

in terms of the parameters in the transfer matrix. Please note that K^ has upper triangular 
form for this choice of parameters. After a displacement of the bosonic operators a — )■ 
a+/3z/]"/2 and a simultaneous rotation of the spin the hamiltonian becomes (up to a constant) 

H = VtQn + A,,S, + g{S^a^ + 5"a) . (4.31) 

Although this operator has been obtained from non-diagonal boundary conditions this oper- 
ator commutes with the charge n — S^. We could not find any choice of parameters leading 
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to a hermitean hamiltonian including counter-rotating terms. The same statement applies 



27 - different from what has been claimed 



to a similar choice of parameters proposed in Ref. 
there. It is worth mentioning that within the above choice of 9 parameters in the transfer 
matrix only four (plus the spectral parameter A) enter the final hamiltonian, which can be 
obtained also from diagonal boundaries. As some of these spurious parameters will be seen 
to enter the TQ-equations, this is a hint that the function Q indeed carries information 
about the eigenstates of the hamiltonian, which are sensitive to local changes of the basis. 
It is furthermore intriguing that despite the presence of a conserved U{1) charge transfor- 
mations such as f l4.12p induced by a transcendental function (see Eq. (I4.16P ) are required in 
order to guarantee solutions to the functional equations that can be parametrized by finitely 
many roots of a polynomial Q{\)- In contrast to the transformation of the factorization of 
the quantum determinant of the full transfer matrix, a suitable t] dependence of the parame- 
ters (see above) is required in the quasi-classical limit in order to ensure that the two lowest 
orders in 7] remain unchanged and that no lower orders are created by that transformation 
in the TQ-equations. Taking account for these subtleties leads to the ansatz 

Q{z) = Hz)Q{z) 

1 \ (4.32) 



:F{z) = ivVT z + 



bJXif 



The quasi-classical expansion of the equation (I4.14p is carried out along the steps described 
in the appendix. The lowest orders l/r/^ and l/r^ of the equation are identically satisfied for 
a; = 1; the order rf leads to the second order differential equation 

^A_,(A)g"(A) + i?(A)Q'(A) + U{\)Q{\) = Ao(A)Q(A) (4.33) 

which has the form of a Schrodinger equation for the Hamilton operator 
n = |A_,(A)92 + R{X)dx + f/(A). 

A_,(A) = -/3^A2(A2-^o') 
R{X) = -^A [2A^ - X\2zl + (3\1 - 2z^ + 2^^)) + P^zl{l + 2z^ - 2^^ 



UiX) ^ -^ 



X\f,-u^l3' + 2{z,-^^-^t)) (434) 

+ ^x' {-AzoU + /5'(3 - 4er + 4^i(i + er)) 

+2^2(2^1 + /ir^r/3' + 2(1 - er - ^t))) - ^^o(i + 4^ier 
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In the present case Aq = Aq^ + aSj^^A^ + A[,^U^. In particular, A[,°^ = U{X = 0), which 
is consistent with a constant solution for Q, as should be expected for a rotating Jaynes- 
Cummings model. Aq is the energy as a function of the spectral parameter; it has to be 
fixed by the requirement that the corresponding differential equation fl4.33p has a polynomial 



solution, in analogy to the parameter ( in Ref. |25|, Eq. (12). 

It is seen that the spurious parameters zi, fi^ , z/f, C,^, and C,q appear in the TQ-equations 
but not in the Hamiltonian. Zi can be safely set to zero, whereas /ij~ and u^ only occur in 
the invariant combination /ij~z/-f A^ = ^i — A^ — ^^^r— ; the same applies to ^^ in terms of 
the trace of K'^. This indicates that their appearance in the TQ-equation reflects a change 
in the eigenbasis encoded in the function Q. 

It is worth noticing that the transformation leading to polynomial Q in the quasi-classical 
limit does not fix the asymptotics of the TQ-equations for the full transfer matrix. This 
highlights that the corresponding models differ considerably. 

V. SUMMARY 

We have performed the separation of variables for spin-boson models generated from a 
bosonic Lax operator descending from rational six-vertex models with twisted and open 
boundary conditions. Our focus was on two-site compositions where a single spin interacts 
with a single bosonic mode. Generic spin-boson interactions are counter-rotating and thus 
do not conserve S^ + n in contrast to the rotating models. This hampers the construction of 
a simple (pseudo) vacuum state, which is the necessary starting point for a diagonalization 
of the model by means of the algebraic Bethe ansatz. We employed the functional Bethe 



Is 



ansatz, proposed by Sklyanin for systems without known reference vacuum state J6|]. The 
application of this method to spin-boson models requires compatibility of the formalism with 
an infinite dimensional bosonic Hilbert space. Our analysis demonstrates that the technique 
of separation of variables - originally designed for finite dimensional representation spaces - 
carries over straight forwardly to this scenario. Specifically, we have found that the infinite 
dimensional bosonic Hilbert space is one-to-one reflected by a half-infinite Sklyanin lattice. 
This provides further indication that the functional Bethe ansatz and its central features 
relying on the Yang-Baxter algebra alone do not depend on the chosen representation. The 
infinite dimensional representation of the 'factorization algebra' A^ [6, Eq.(3.13)) is due to 
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the absence of a bosonic highest weight and leads to a pecuhar factorization of the quantum 
determinant, which in turn is determined to some extent by the boundary d& = {x~, xf}x 
{x^} of the Sklyanin lattice &. As for finite dimensional representations (e.g. spin models), 
the Sklyanin lattice results to be composed of eigenvalues of the spin and boson part of the 



operator zeros of an off-diagonal monodromy matrix element 49|. Interestingly enough, a 
suitable gauge transformation of the monodromy matrix has facilitated the diagonalization 
of these operator zeros for the twisted model. 

The TQ-equations appear as a Schrodinger equation for the action of the transfer ma- 
trix on the Sklyanin lattice, where Q is formally related to the eigenstates of the model 
hamiltonian by some isomorphism [6|. Algebraic Bethe equations are obtained from the 
TQ-equations for polynomial Q. However, the function Q is by no means uniquely defined. 
A TQ-equation with a transformed function Q can be obtained in general after transferring 
the effect of a transcendental part of Q onto the factorization of the quantum determinant. 
In terms of the non-hermitean operator zeros x this induces a renormalization of their right 
(and left) eigenstates, which in turn is mathematically equivalent to the isomorphisms for the 
factorization algebra A'a described in 6|. We have analyzed this ambiguity for the function 
Q on a class of spin-boson hamiltonians with the aim to arrive at a TQ-equation allowing 



for a polynomial solution Q. We find that the 'canonical' factorization [50| of the quantum 
determinant of the transfer matrix will typically deviate from the proper factorization 51 1 
even in this simple case. 

A presumable physical meaning of this proper normalization remains unclear and would 
be worth a further investigation. There are two pieces of evidence for a physical impact 
behind this choice of normalization. One is that different transformations are needed for 
the full transfer matrix and its quasi-classical limit. The second is that the Bethe equations 
are considerably changed in a way shown in Eq. (I4.18p . Understanding these issues would 
constitute a significant step forward in the theory of integrable quantum systems. 
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Appendix A: Quasi-classical limit 

The so called quasi-classical limit of the transfer matrix |7|, |35| consists in a series ex- 
pansion in the 'quantum parameter' rj (playing the role of h) of the transfer matrix around 
r] = 0: i{\) = f*^°^ + 7] f^-^'^ (A) + r/^ f*^^^ (A) + . . . with the aim of creating a commuting family 
of quasi-classical transfer matrices r*^^)(A). This procedure has proved to be particularly 
useful for extracting 'simple' though non-locally interacting hamiltonians out of the transfer 



matrix. Examples are the Gaudin magnets and corresponding BCS-like models. 36 



-|39| 



There is much freedom for introducing an ?7-dependence to the boundary matrix parame- 



ters in integrable theories. Examples for the quasi-classical 
and open boundary conditions can be found in Refs. |4Ch42 



imit of spin models with twisted 



In the present cases (for twisted or open boundaries) the transfer matrix is a finite sum: 

t(A) = r^-'^f^-"') + ■ ■ ■ + f(°) + r] f(^) (A) + r]^ f^^) (A) + ■ ■ ■ + r/™ f^™) (A) (Al) 

where k and m are integers. Expanding the commutator relation [t(A),t(A')] = in 77, we 
obtain 

2m 

[t(A),t(A')]= 5^ VQ(A,A') = 0, 

l=-2k 

which implies C;(A, A') = for all /. The first relevant terms are 

C_2fc+i(A,A') =[f^''\\),f^-'^'\\')] + [f^-'^'\\),f(-'\\')] , 

C_2fcH-2(A, A') =[f(~'=)(A),f(-^+2)(y)] ^ [f(->^+^\X),T^->'\\')] + [f(-'=+i)(A),f(-'=+i)(A')] , 

(A2) 

From the expressions above, one finds that the first f'^")(A) which is not a C-number (times 
the identity) gives rise to a family of commuting operators. Generically, the lowest or- 
der f'-"^-' is a C-number. Therefore the first class of integrable models are generated by 
[f(-'=+i)(^_\), f'-~'^"'"^''(A')] = 0. In the presence of boundary matrices, these are typically non- 
trivial operators but representing non-interacting hamiltonians. The task is then to tune 
the free parameters such that f^^'^^^^X) is also a C-number. The lowest non-trivial order 
in T], e.g. f*^^^^^^(A), is typically a hamiltonian with non-trivial interactions. 
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The TQ-equation (14.91) is a second order difference equation which has to be solved to 
determine the spectrum of the spin-boson system. In the quasi-classical limit 77 — )■ it 
becomes a differential equation |7[ for the Q-function we want to study here. 



m ra / 2 \ 

j=-k j=-k ^ ^ 

m , 2 \ 



(A3) 



j=-k 

where f ^ are suitable factorization of the quantum determinant (defined above as 'B.'^IS^ or 
A''^ for quasi-periodic or open boundaries respectively). We point out that the differential 
equation arising from the n-th order of the ?7-expansion involves (solely) the coefficient A„ 
of the eigenvalue; therefore, without the knowledge of the latter, the differential equations 
arising from the TQ-equation can be integrated only at a formal level. 
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